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ABSTRACT 

The Suzaku observation of LMC X-3 gives the best data to date on the shape of the accretion 
disk spectrum. This is due to the combination of very low absorbing column density along this 
line of sight which allows the shape of the disk emisison to be constrained at low energies by the 
CCD's, while the tail can be simultaneously determined up to 30 keV by the high energy detectors. 
These data clearly demonstrate that the observed disk spectrum is broader than a simple 'sum of 
blackbodies', and relativistic smearing of the emission is strongly required. However, the intrinsic 
emission should be more complex than a (color-corrected) sum of blackbodies as it should also 
contain photo-electric absorption edges from the partially ionised disk photosphere. These are 
broadened by the relativistic smearing, but the models predict ~ 3-5 per cent deviations for 1/3- 
1 solar abundance around the edge energies, significantly stronger than observed. This indicate 
that the models need to include more physical processes such as self-irradiation, bound-bound 
(line) absorption, and/or emission from recombination continuua and/or lines. Alternatively, 
if none of these match the data, it may instead require that the accretion disk density and/or 
emissivity profile with height is different to that assumed. Thus these data demonstrate the 
feasibility of observational tests of our fundamental understanding of the vertical structure of 
accretion disks. 
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It is important to study the emission from 
the accretion disk to understand how the gravita- 
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tional energy is converted to radiation. The sim- 
plest models of accretion disk spectra assume that 
the gravitational energy dissipated at each radius 
thermalizes to a blackbody spectrum. Summing it 
over all radii, under an appropriate inner bound- 
ary condition, produces the well-known Shakura 
& Sunyaev (1973) disk model. This is derived in 
Newtonian gravit y, but was extended to fu ll gen- 
eral relativity by IHovikovAThomi (jlili)- The 
maximum temperature produced by a disk accret- 
ing close to the Eddington limit around a 10 Mq 
black hole is 1-2 keV, easily observable with X- 
ray satellites, whereas a similarly Eddington lim- 
ited Active Galactic Nuclei (AGN) of lO'^M© has 
a temperature of 30-60 eV, in the unobservable 
EUV regime. Thus stellar remnant black hole bi- 
nary (hereafter BHB) systems give a much better 
test of disk models than AGN. 
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The emission from each radius is a true black- 
body only when the disk is effectively optically 
thick to absorption at all frequencies. Free-free 
(continuum) absorption drops as the frequency in- 
creases, so the highest energy photons from each 
radii are unlikely to thermalize. This forms in- 
stead a modified (or diluted) blackbody, with ef- 
fective temperature higher than that for complete 
thermalization by a factor of /coi (termed a color 
temperature correction) . The full disk spectrum is 
then a sum of these modified blackbodies, but this 
can likewise be approximately described by a sin- 
gle color temperature correction to a sum of black - 
body disk spectrum dShimura fc Takaharal llQQSl ). 
However, continuum processes do not necessarily 
dominate the total absorption at all frequencies. 
Bound-free (photo-electric) absorption from par- 
tially ionised metals can be important, especially 
at high frequencies where th e free-free absorp - 



tion becomes less significant ( Davis et al. l2005l) 



This imprints atomic features onto the emission 
from each radius, distorting the spectrum from a 
smooth continuum. The strength of these features 
is set by radiative transfer through the vertical 
structure of the photosphere, so they are one of 
the few diagnostics o f the internal disk properties 
(|Done fc Da.visll2008l ). 

Thus the intrinsic spectrum from each radius 
can be complex, but each of these is smeared out 
by the combination of special and general rela- 
tivisic effects which arise from the rapid rotation 
of the emitting material in a strong gravitational 
field (jCun ningham f975). The resultant smeared 
spectra are summed together to form the total disk 
emission, which is not that different from a sum 
of smea red color ternperature corrected blackbody 



spectradDavis et al 



20051 iDavis fc Hubenvl [20oi 



Done fc David 2008 ). We need excellent data in 
order to detect these smeared spectral features, 
and use them to test our understanding of the disk 
vertical structure. 

The BHB LMC X-3 is the best object currently 
known for this. It has the lowest absorption col- 
umn density of any BHB which shows disk domi- 
nated spectra, extending the bandpass over which 
the data can constrain the models down to the 
softest X-ray energies. Additionally, the system 
parameters are fairly well determined, especially 
distance which is known to better than 10 per cent 
due to its location in the Large Magellanic Cloud. 



Thus the conversion of flux to luminosity is subject 
to smaller uncertainties than for any other BHB. 
This combination of properties means LMC X-3 
offers a unique laboratory for testing our under- 
standing of accretion disk physics. 



Here we use data from a Suzaku ( Mitsuda et al 



20071 ) observation of the disk dominated state 
in this object. This satellite combines moderate 
spectral resolution CCD's useful to constrain 0.7- 
10 keV spectrum of the source, together with the 
Si PIN photo-diodes of the non-imaging Hard X- 
ray Detector (HXD, Takahashi et al. 2007) which 
can simultaneously determines the 12-30 keV spec- 
trum. Thus these data give a complete picture of 
the emission, being able to constrain the (weak) 
tail to higher energies which otherwise gives a sig- 
nificant source of uncertainty in reconstructing the 
disk spectrum from CCD data alone. 

The spectrum of LMC X-3 during this observa- 
tion is dominated by a clear disk component with 
luminosity around 10 per cent of the Eddington 
limit. This is low enough to avoid the multiple 
uncertainties which arise at higher luminosities, 
where the disk may puff up due to radiation pres- 
sure, advection of radiation may become impor- 
tant, and strong winds from the inner disk can 
distort the spectrum both through absorption and 
from changing the intrinsic disk spectrum due to 
the mass loss. Thus these observations give the 
most sensitive test to date of our understanding of 
a 'clean' accretion disk to compare with the mod- 
els. 

2. Observation and Data reduction 

We observed LMC X-3 with Suzaku from 
2008 December 22 UT07:14 through December 23 
UT20:49 (epoch 5). Suzaku carries 4 X-ray tele- 
scopes (XRT, Serlemitsos et al. 2007), each with 
a focal-plane X-ray CCD camera (X-ray Imaging 
Spectrometer; XIS, Koyama et al. 2007) operat- 
ing in the energy range of 0.2-12 keV. Three of 
the XIS (XISO, 2, 3) have front-illuminated (FI) 
CCDs, while XISl utilizes a back-illuminated (BI) 
CCD, achieving an improved soft X-ray response 
but poorer hard X-ray sensitivity. Since XIS 2 is 
no longer available, we use the two remaining FI- 
CCDs (XIS03) together with the BLCCD (XISl). 
We combine these with data from the Si PIN 
photo-diodes of the HXD (|Takahashi et al.ll2007t ). 
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which covers the 10-70 keV energy band. 

The source was centered in the XIS field of view. 
The XIS were operated with the 1/4 window op- 
tion so that pileup did not become significant. We 
use version 2.2.11.22 of the pipeline provided by 
the Suzaku team, and the XIS data are screened 
using the standard criteria (i.e. only GRADED, 2, 
3, 4, 6 events are accumulated, time interval af- 
ter passage through the South Atlantic Anomaly 
is larger than 436 seconds, the object is at least 5° 
and 20° above the rim of the Earth during night 
and day, respectively). This gives a net exposure 
of 73.98 ks. 

We extract XIS events from both 5x5 and 3x3 
editmodes using a circular region with a radius of 
4'. 3 centered on the image peak. This is larger 
than the window size so the effective extraction re- 
gion is the intersection of the rectangular window 
with this circle. We use ftool xisrmfgen (ver- 
sion 2007-05-14) and xissimarfgen (version 
2009-0 l-08)to calculate the response matrix and 
auxiliary response of the instrument. We co-add 
the data from the two FI CCDs as these are very 
similar, and generate the co-added response for 
this XIS03 spectrum using the ftools ADDRMF 
and ADDARF. We add systematic errors of 1% to 
each energy bin for both XIS03 and XISl spectra. 

The source is bright, with count rate from XISl 
of 24 cts s^^ in the range of 0.7-10 keV, which 
corresponds to 3.1 x 10^^" erg s~^ cm~^. Though 
background can be important towards the highest 
energy end of this bandpass, the large point spread 
function of Suzaku means that the 1/4 window 
data cannot be used to simultaneously determine 
the background as there is no source-free region. 
Instead we use the Lockman hole data obtained by 
Suzaku from 2009 June 12 UT 07:17 through June 
14 UT01:31 with an exposure of a net exposure of 
92 ks, using the same shaped (intersection of circle 
and rectangle) region. The background levels are 
negligible in soft energy band, and reach to ~0.7 % 
and - 5 % of our LMC X-3 data at 8 keV, for 
XIS03 and XISl, respectively. In this paper, we 
used the 0.7-10 keV band and the 0.7-8 keV band 
for XIS03 and XISl, respectively. We show the 
background subtracted XIS03 and XISl spectra 
in Fig. m 

The HXD was operated in the nominal mode 
throughout the observations. We use the standard 
pipeline processed PIN data (object is at least 



5° above the earth rim, time interval after pas- 
sage through the South Atlantic Anomaly is longer 
than 500 seconds, and cutoff rigidity is greater 
than 6 GV). We extract background from the 
same goodtime intervals of the PIN background 
model based on lcfitdt method (2.0ver0804) 
provi ded by the HXD tea m for each observa- 
tion ( Fukazawa et al.l [20091 ). Here we note that 
the HXD-PIN background model is aimed to re- 
produce the non X-ray background. The source 
exposure is corrected for the 7.3 % deadtime, 
which is calculated from the pseudo event r ates 



(|Takahashi et all 120071 : iKokubun et"al1 l2007l ) by 



using the ftool hxddtcor, giving a net ex- 
posure of 70.03 ks. This is ^ 25% less than 
the XIS, but the source shows little variability 
so we do not exclude any XIS data in order 
to have only exactly simultaneous data. These 
data show that the source is not significantly de- 
tected at the highest energies, so we restrict the fit 
range to 12-30 keV using the appropriate response 
ae_hxd_pinxinome5_20080716.rsp. 

The background subtracted PIN spectrum is 
shown in Fig. [T] together with the XIS spectra. 
This has a 12-30 keV count rate of 0.076 cts s^^. 
This corresponds to a flux of 2. 2x 10~^^ erg s^^ cm" 
which is only ^ 4 times higher than flux of the 
cosmic X-ray background (CXB) in the same en- 
ergy range. Thus the CXB is not negligible, and 
we include t his as a fixed exponentially cutoff PL 
(|Boldtlll987| ) in all our spectral modeling for the 
PIN dat43- 



-^http: / /heasarc. gsfc.nasa.gov/docs/suzaku/analysis/pin_cxb. html 




Energy (keV) 



Fig. 1. — The background subtracted X-ray spec- 
tra of LMC X-3 obtained with XIS03 (black), 
XISl (red) and the HXD PIN (green). 
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3. Analyses and Results 

We describe here the system parameters of 
LMC X-3 adopted in this paper. LMC X- 
3 is a persistent source, so the emission from 
the outer accretion disk always contributes to 
the optical spectrum, and X-ray irradiation 
can alter the structure and spectrum of the 
companion star. This introduces some un- 
certainty in determining the properties of the 
comp anion, and hence the mass of the black 
hole (van der Klis et al. 1985 : Cowley et al. 1983t 



Soria et al.ll200ll: iNegueruela fc Coell2002[) . Tak- 
ing account the uncertainties, we consider the 
black hole mass in the range of 7-9Mq and incli- 
nations of 50° -67°. Thus for each physical model 
(kerrbb and bhspecH in xspec) we report fit 
values over this range, showing also results for the 
specfic combination of 7AIq and inclination of 67° 
used bv iDavis et al. I 12006). We assume a source 
distance of 52 kpc ( di Benedettol 1997 ). 

We use the xspec spectral fitting package. 
Large fit residuals due to calibration uncertain- 
ties are often observed near the edge structures 
of the XIS/XRT instrumental responses. We ex- 
clude energy range of 1.6-2.0 keV and 2.2-2.4 keV, 
and model the remaining instrumental feature at 
~ 3.2 keV by gold M edge by a gaussian line 
wi th fixed gauss i an w idth of 0.1 keV as is done 



Kubota et al.l ( 20071) 



The absorption column de nsity is well de- 
termined as 3.8 x 10^° cm-^ (|Page et al.l l2003f ) 



from the depths of edges measured by the XMM- 
Newton Reflection Grating Spectrometer data. 
This is consistent with most of the absorption 
being due to our own Galaxy as the colum den- 
sity along this line of sight is 3.2 x 10^° cm~^ 
asmeasured by radio observations ( Nowak et al.l 
20011). However, low-energy calibration (especially 
the quantum efficiency) of the Suzaku XIS suf- 
fers from relatively large systematic errors due to 
the uncertainties of the contamination thickness 
and its composition. To minimize the impact of 
these systematic errors while keeping good low- 
energy coverage, we restricted the energy range 
above 0.7 keV in the spectral analysis. In this en- 
ergy range with the calibration data base (CALDB: 
VERSION XIS20090203), we estimate the system- 



atic error in determining TVjj is about a few times 
10^" cm~^ (XIS team, private communication). 
Because this is comparable to the column den- 
sity to LMC X-3, we allow (described using 
PHABS) to be free in the spectral analysis unless 
otherwise stated. 

When we analyze the disk emission, it is im- 
portant to model the power-law component ap- 
propriately, because its low-energy part sometimes 
affect the fitting of the disk emission. A simple 
power law (hereafter PL) plus CXB fit to the PIN 
data alone in the 12-30 keV band gives a pho- 
ton spectral index of 1.66-2.38 with the best fit 
value of 2.01. Because the power-law component 
in the high/soft state usually has a photon index 
larger than 2.0 (e.g.. Done et al. 2007 and refer- 
ences therein), we consider that the true index of 
LMC X-3 lies between 2.00-2.38. Therefore, we 
constrain the index in this range in the spectral 
analysis rather than allowing it to be free. Oth- 
erwise, the power-law component would be opti- 
mized to fill the fit residuals in the low-energy 
band, where the statistics are much better than 
those of the HXD PIN, and the best-fit power- 
law index could significantly deviate from the true 
slope in the hard energy band. 

3.1. Empirical modeling of the disk: DISKBB 

In order to compare the spectral data with the 
previous observations, we first use the common 
black hole spectral model of an absorbed disk plus 
power law. The simplest disk model is DiSKBB, 
which has a temperature distribution T(r) oc 
r~^/ ^, i.e. has no stress free inner boundary cond i- 
tion (jMitsuda et al.lll984t iMakishima et al.lll986l) . 



Figure [2] shows the XIS and the HXD/PIN 
spectra with this best fit diskbb plus PL model, 
with parameters given in table [TJ The spec- 
trum is well fit ix^/dof = 770.25/706) by 
the dominant DISKBB component with Tin of 
0.84 keV. The 0.7-30 keV flux is estimated as 
3.30 X 10~^° erg cm~^, which gives an ab- 
sorbed X-ray luminosity as 1.07 x 10"^^ erg for 
an isotropic emission at D — 52 kpc. This lumi- 
nosity is about 10% of the Eddington luminosity 
for a 7 Mq black hole, and it is the same level 
as that seen in the dimmer of the two ASCA ob- 



^http: / /www. sns.ias.cdu/'^swd/xspec. html 



servations reported by (jKubota et al.l 120051 ) . The 
radius derived from the DiSKBB normalization, 
r?^cosi/{D/10 kpc)2 = 26.0 ± 0.3 is also consis- 
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Table 1 

The best fit parameters based on the diskbb model 



model 


TVh 
10™ 


keV 


disk norm 


r " 


norm or / ^ 




diskbb + PL 
simpl*diskbb 


6.04 ± 0.02 
< 7 X 10"^ 


0.840 ± 0.002 

u. '»'_0.001 


26.0 ± 0.3 

42 8+"-^ 


2.380tp 
2.380_(, 003 


U.UCiO't_o 0006 
„ , , I.+0.0007 

^'■-'--^^■'-n.0023 


770.25/706 
1336.14/706 



Note. — Errors represent 90 % confidence. 

^normalization is defined as r'^^ cos i / {D / IQ kpc)^. 

^Value of plioton index F is limited to be 2.0—2.38 (see text). 

^PL normalization at 1 keV in the unit of photons s^^ keV^^ or scattering fraction / of the SIMPL model. 




Fig. 2. — vFy spectrum of LMC X-3 from the best 
fit DISKBBB plus PL model (a), residuals between 
the data and DISKBB plus PL fit (b) residuals be- 
tween the data and best fit SIMPL*diskbb model. 
In panel (a), the CXB component is subtracted 
from the PIN spectrum. 



tent with the values of 25 ± 3 and 24 ±2 obtaine d 
from the ASCA observations ()Kubota et al.ll2005[) . 

Even though the fit is good, a closer inspection 
of the figure shows that this is not a good physi- 
cal description of a model where the disk provides 
the seed photons for Compton upscattering into 
the power law tail as the PL extends below the 
disk at low energies. Hence we replace the PL 
model with the SIMPL (jSteiner et alj 120081) model 
for Compton upscattering. This takes some frac- 
tion of the disk seed photons and upscatters these 
to higher energies. It is a convolution model, so 
requires an extended energy range which we fix at 
0.1 keV to 1000 keV. 

The fit becomes substantially worse with this 
more physical model {x^/dof — 1336/706), as 
shown by the residuals between the data and the 
best fit SIMPL*DISKBB model (Fig. In ad- 

dition, A'h is much smaller than expected, even 
though it is comparable to the current systematic 
uncertainty of the A'h determination with the XIS. 
Both these effects show that the observed disk 
shape gives the excess at lower energies. In other 
words, a peak profile of the observed disk spectrum 
corresponding to the maximum disk temperature 
is much broader than predicted by diskbb, mo- 
tivating us to do more detailed spect ral analyses 
Simila r conclusions were reached by iDavis et al 



l)2006h . who analyzed BeppoSAX data. Although 
their data had poorer signal-to-noise, it had simi- 
lar low energy coverage and they also found poor 
fits with DISKBB once the PL was reasonably con- 
strained. 
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Table 2 

The best fit parameters based on the kerrbb model 



mass 
Mo 


i 

degree 


10^° cm-" 




M '' 
10'** g/s 


r 


/ 


xVdof 


7 
7 
9 
9 


67 
50 
67 
50 


0.16 ± 0.14 
0.8 ± 0.1 

0.31-,; 
1.0 ± 0.2 


-0.07 ± 0.01 
n +0.010 

^■^^-0.009 

0.738 ± 0.007 


3.32 ± 0.03 
1.27 ± 0.01 
2.58 ± 0.02 
1.005 ± 0.009 


2.38lo Q4 

„ ,7+0.01 
^■■"-0.10 

2-38_o 06 

„ o, +0.07 

^"'-^-o.oo 


n 093+" '"'" 
u.uao_o 004 

U.088_o 008 
U.U»l_o 005 

U.082_o 007 


783.58/706 
772.86/706 
776.33/706 
776.70/706 



Note. — Same as table ^ but for the KERRBB model. 

'^Blackhole spin parameter. 

^Mass accretion rate of the KERRBB model. 



3.2. EfTect of the general relativity: tem- 
perature profile and relativistic smear- 
ing : DISKPN and kerrbb 

We now use progressively more complex mod- 
els for the disk emission together with the SIMPL 
Compton upscattering model for the tail. First, 
we consider the diskpn model which includes 
the stress-free inner boundary condition assuming 
a Paczynski-Wiita (pseudo-Newtonian) potential. 
This smoothly connects the temperature profile 
from T(r) cx r"'^/'* at large radii, and T(rin) = at 
Tin, where rjn = GRg, in a way which is similar to 
that expected from the ful ly relativistic equation s 



for a zero spin black hole ( Gierlinski et al. 19991 ) 



This has the effect of broadening the spectrum 
slightly as it reduces the emission from the hottest 
material, and hence gives a better fit to the data 
than the SIMPL*diskbb {xl = 984/706), though 
still not as good as the unphysical DiSKBB-fPL 
model. This indicates that the soft component in 
the data is still broader than the diskpn descrip- 
tion. Relativistic corrections to the temperature 
profile alone do not sufRciently broaden the spec- 
trum to match the observed data. 

The next step up in complexity is to include the 
special and general relativistic effects which dis- 
tort the observed sp ectrum at each radius. We use 
the KERRBB model ( Li et ahlboOSl) for this, which 
calculates the intrinsic temperature distribution 
from a fully relativistic, stress-free inner boundary 
condition for arbitrary spin, multiplied by a color 
temperature correction factor, /coi, to approxi- 
mate the effect of modified blackbody emission. 
We fixed /cd = 1.7 as appropriate for the XIS 



CCD bandpass at this luminosity (Done & Davis 
2008). This intrinsic emission from each radius is 
convolved with the relativistic transfer function for 



that radius and assumed disk inclination! Li et al 



20051). 



We explore the dependence on mass and incli- 
nation by fitting the model with different combi- 
nations of these two parameters. These fit results 
are shown in table [2] and Fig. |3l These all give 
a fit which is similarly good as that derived from 
the DISKBB-I-PL model (see table [5]). The quality 
of the fit is not very dependent on mass and incli- 
nation, but the derived value of the spin changes 
from almost zero (lowest black hole mass, highest 
inclination) to ~ 0.7 (highest black hole mass and 
lowest inclination). 

We illustrate the fit in Fig. [3^, showing the 
vFi, spectrum and residuals for the models with 
mass of 7M(7) and inclination of 67°, as used by 
iDavis et all ( 2006 ). This gives a good fit {xl = 
776/706) with flat residuals across the whole en- 
ergy band (see Fig. I^b), and gives a spin of a* ~ 
0, con sistent with that derived by iDavis et~al 
( 20061 ). Relativistic effects on the disk contin- 



uum were first proposed bv ICunninghaml (1975), 
but this is the first time they have been signifi- 
cantly detected, as it is much harder to disentan- 
gle this smearing on a bro ad continuum than on 



a line ( Fabian et al. 19891 ). Though the smaller 



values of iVn may suggest that the data still has 
excess soft emission and is somewhat broader even 
than that produced by the temperature profile and 
relativistic broadening, the calibration uncertain- 
ties do not enable us to conclude unambiguously 
whether there are other broadening factors. 
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3.3. Effect of the vertical structure of the 
disk : BHSPEC 

The assumption of a constant /coi is only an 
approximation to the true spectrum from a given 
radius. In the real disk, changes in opacity at 
atomic edges produce spectra intrinsicahy broader 
than a blackbody. Under the assum ption of the 
solar metallicity, the bhspec model (jPavis et al 



20051 ) uses radiative transfer to calculate the in- 



trinsic spectrum at a given radius, so does self- 
consistently includes the color temperature cor- 
rection. It then convolves this with the relativistic 
transfer functions to produce the best disk spectra 
to date. 

We fit the spectra with the bhspec disk model 
for the same combinations of i and M as for the 
KERRBB fits above. The results are shown in ta- 
ble [3l The model gives similar spins to those de- 
rived from KERRBB (see tabled]), showing that the 
assumption of /coi = 1.7 is appropriate to the 
KERRBB fit for the XIS CCD bandpass. More- 
over, the estimated absorption is now consistent 
with the known column density in this direction ^ 
3 X 10^° cm~^. Thus the bhspec model produces 
additional excess soft emission than produced by 
relativistic broadening alone. However, despite 
being a more physical model, it gives a much worse 
fits than KERRBB (x'^/dof = 1212/706 for the case 
of M = 7Mq and i = 67° compared to 783/706). 
Figure 13] shows the residuals between the data and 
the model for each mass and inclination. These 
are all clearly dominated by a broad feature be- 
low 1 keV, which is caused by the smeared atomic 
absorption edges present in the model. 

We checked the result by including relativis- 
tically smeared ionised reflection, modelled us- 

ing the convolution v ersion of the models of 

iBallantvne et al. ( 2001 ) as described in Done fc Gierlihsldl 
( 20061 ) ■ Though the reflection itself cannot be well 
constrained because of couplings to mass, incli- 
nation and the difference of the continuum disk 
model, it does not change this conclusion: bhspec 
gave a worse fit to the data than kerrbb if we 
did not constrain its absorption, and the bhspec 
residuals are dominated by features around 1 keV. 




Energy (keV) 

Fig. 3. — vFi, spectrum of LMC X-3 based on the 
best fit SIMPL*kerrbb model (a), and residuals 
between the data and the model (b). Same as 
Fig. 12] the CXB subtracted spectrum is shown for 
the PIN data (panel a). 



(a) M=7Mo, i=67' 



|I^Tl''':,ll|ii II Wf i 




5 10 20 

Energy (keV) 

Fig. 4. — Residuals between the data and the best 
fit bhspec models with different mass and incli- 
nation. (a)M = 7Mq and i = 67°, (b)M = 7Mq 
and i = 50°, (c)M = 9Mq and i = 67°, and 
(d)M = 9Mq and z = 50°. 
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Table 3 

The best fit parameters based on the bhspec model 



mass 

Mq 


i 

degree 


10^° em-^ 


a* 


10^** erg/s 


r 


/ 


xVdof 


the solar abundance 


7 
7 

9 


67 
50 
67 


2.7 ± 0.1 
3.9 ± 0.1 
3.1 ± 0.1 


0.542 ± 0.006 
0.286 ± 0.009 


1.799 ± 0.005 
1 054+""''^ 
1.748 ± 0.005 


2.380;!;q 

2.380_Q QQg 
2.380_Q QQg 


102+"™^ 

U.lUZ_o 002 

0.096 ± 0.002 
0.100 ± 0.001 


1212.47/706 
1058.56/706 
1094.30/706 


9 


50 


4.3 ± 0.1 


0.742 ± 0.004 


1 085+"°°^ 

l.UOO_o 


2.380_(, 


0.094 ± 0.002 


974.65/706 


1/3 solar abundance 


7 
7 


67 
50 


1.6 ± 0.1 

2.7 ± 0.1 


-0.076 ± 0.010 
0.536 ± 0.006 


1.033 ± 0.002 


2.380+0 007 
2.380_o.oo8 


0.102 ± 0.002 
0.098 ± 0.002 


1151.91/706 
1062.58/706 


9 


67 


2.0 ± 0.1 


0.278 ± 0.009 


1.717 ± 0.005 


2.380_(, 006 


0.102 ± 0.002 


1063.26/706 


9 


50 


3.1 ± 0.1 


0.737 ± 0.004 


1.061 ± 0.002 


2.380_Q 009 


0.096 ± 0.002 


987.75/706 


1/10 solar abundance (just for comparison) 


7 
7 
9 
9 


67 
50 
67 
50 


0.7 ± 0.1 
1.69 ± 0.1 

2.0 ± 0.1 


-0.049 ± 0.010 

+ 0.007 

0.748 ± 0.004 


. -.„_+U.lJlJl3 

'^''-0.004 

1 011+"'"'^ 
1.680 ± 0.005 
1 039+""''^ 

J-.uo»_Q QQ3 


2.380;';o 008 
2.380_o 008 
2.380_Q 009 
2.380^0 010 


0.102 ± 0.002 
0.099 ± 0.002 
0.101 ± 0.002 
0.097 ± 0.002 


1014.46/706 
974.17/706 
973.39/706 
919.16/706 



Note. — Same as tableland tablc[2]but for the BHSPEC model. 



3.4. Constraint on the absorption lines 
due to highly ionized iron 

Absorption lines from highly ionized Hydrogen 
(H-) and Helium(He-) like iron are often seen in 
disk dominate d black hole bina ries (see e.g. the 
compilation in iDone et alJboOTt) . We focus on the 
iron K line region to sensitively search for such 
features. We fit the 5-8 keV data with a power 
law and narrow absorption line (with width fixed 
at 10 eV). There is evidence for such a line at 
6.61±0.04 keV, of equivalent width of — 15±6 eV. 

Though this is significant at about the 99% con- 
fidence level on an F-test as xt drops from 233/221 
to 215/219, the line energy is slightly too low to 
be He- like iron at 6.70 keV, even if we consider the 
calibration uncertainty of the gain, ^ 20 eV, when 
the 1/4 window option is usec^. The difference 
may be even larger, because the lines are often blue 
shifted due to the disk wind. As seen i n the Chan- 
dra/H ETGS spectra of GX 13-H (fUeda et al 

' 20091) . He- 



200J and GRS 1915-^105 ( Ueda et al. 



like iron energy can be observed at low er energies 
than the prediction. Ueda et al. ( 20091 ) suggested 
that the slight shift of the line energy can be 
caused by a possible contamination of less ionized 
ions such as Li-like ones at 6.68 keV, uncetainties 



^Based on the CALDB version XIS20090203 



in the incident line energies for He-like (and more 
electrons) iron, and/or the velocity-field structure 
in the wind along the line-of-sight. However, even 
considering these effects and the XIS calibration 
uncertainty, 6.61 ± 0.04 keV is still difficult to be 
interpreted as He-like iron absorption lines. Fixing 
the line energy at 6.70 keV gives xt = 226/220, so 
again formally significant at the ~ 99% level, with 
an equivalent width of —9 ± 6 eV. This value rep- 
resents the conservative constraint on the He-like 
iron absorption line. 

Other BHB at similar continuum luminosi- 
ties 10^^ erg s~^) show absorption lines of 
equivalent width of ^ 30 eV in b oth H- and 
He-like iron ( e.g., GRO J1655 - 40 (lUeda et al 



1998t lYamaoka et"all l200lt iMiller et al.l |2006[) . 



IH 1743 - 322 (|Miller et al.ir2006bf )1. These BHBs 
all have similar inclination angle of i ~ 70°. Thus 
less significant absorption lines may indicate that 
the inclination is likely to be lower than 67°. Al- 
ternatively, these BHBs showed higher disk tem- 
perature Tin than our LMC X-3 data, typically 
at > 1.0 keV when they showed absorption lines, 
and thus the disk temperature may also be impor- 
tant in creating the disk wind in addition to the 
luminosity. 
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4. Discussion & Conclusions 
4.1. Black hole spin 

The resulting values of spin are fairly similar for 
both BHSPEC and kerrbb with /coi = 1.7, giving 
a range in spin from 0-0.7 depending on mass and 
inclination. Although this seems a large range, it 
only corresponds to a change in inner radius by 
less than a factor of 2, from 6Rg at a* = to 
3ARg at a* = 0.7, whereas for extreme spin at 
a* = 0.998 the inner radius becomes as small as 
1.23Rg. Thus the models are not very sensitive to 
spin for low to moderate spin, even for the best 
constrained black hole system parameters, as the 
implied change in radius for the last stable orbit 
is low. However, they are very sensitive to high 
spin, and strongly rule out such values for LMC 
X-3. 

Hence we conclude that the black hole in 
LMC X-3 has low-to-moderate spin. This 
is also the case for all other spin determina- 
tions from disk continuum fits (Davis et al. 2006t 



Shafee et al]l2006t IGou et al.ll2Q09l ). apart ^rhaps 
from GRS 1915 + 105(McClintock et al 2006 but 
see Middleton et al 2006). Such low to moder- 
ate spins are in line with current (but probably 
quite uncertain) theoretical models which indi- 
cate birth spins of 5, 0.8 — 0.9 from stellar collapse 
( Gammie et al ] |2004|) . Jowever, we note that the 
iron line fits often imply rather higher spins, some- 
time for the same objects (e.g. Miller et al 2009). 
This is clearly an area of cu rrent intense research, 
see e. g. the discussion in iKolehmainen fc Pond 
(|2009[ ). 



4.2. Vertical Structure 

The Suzaku observation of LMC X-3 shows 
a spectrum which is dominated by the accre- 
tion disk, giving the best data yet on its de- 
tailed spectrum, and a simultaneous determina- 
tion of the high energy tail. We model the tail as 
Compton upscattering of seed photons from the 
disk. The disk spectrum is much broader than 
a simple sum of blackbodies with T(r) cx r^"^/^ 
(diskbb). It is also broader than a sum of black- 
bodies with temperature profile given by the (ap- 
proximate) stress-free inner boundary condition in 
the pseudo-Newtonian potential (diskpn). Rela- 
tivistic smearing of the intrinsic blackbody spectra 



from each radius (with temperature profile given 
by the stress-free inner boundary condition) is re- 
quired before the model gives a good description 
of the overall shape of the data (kerrbb), though 
the kerrbb model also requires an absorption col- 
umn close to zero. 

The BHSPEC model dispenses with the color 
temperature correction by calculating the spec- 
trum self-consistently from radiative transfer 
through the vertical structure of the disk. The 
column density derived from this model is consis- 
tent with the known column, but the fit is poor 
due to residuals at 1 keV. In order to clarify the 
reason for the residuals seen in the bhspec fit, 
Fig. [S] shows a comparison of the best fit bhspec 
model with that of kerrbb. The bhspec model 
gives more soft emission below 1 keV, and it does 
this due to the presence of absorption edges, most 
importantly the K edges of H-like Nitrogen and 
Oxygen at 0.66 keV and 0.87 keV, respectively. 
Though these are broadened and smeared by the 
relativistic effects, they are still deeper than any 
features present in the data, so appear in the resid- 
uals as an excess at the absorption edge energy. 
These residuals are around 5% (peak to peak), but 
the data are actually consistent with no structure 
at these energies, as seen by the featureless resid- 
uals from the best fit KERRBB model (Fig. |3b). 

To reduce the size of the atomic features, we 
first consider sub-solar abundances, since the 
LMC i s suggested to have 1/2 ^1/3 solar abun- 
danc e ( Russell fc Bessell 1989t Russell fc Dopita 
1990l : Korn et al. 2002[ ). In order to examine this 
effect, we extended the original BHSPEC model 
to account for sub-solar metallicity. We fit the 
data with the new bhspec disk model of fixed 
1/3-solar metallicity for the same combinations 
of i and M as used previously. The results and 
fit residuals are shown in table [3] and Fig. El re- 
spectively. For the condition of M = 7 Mq and 
i = 67°, the best fit bhspec model with 1/3- 
solar abundance is shown in Fig. [5] together with 
the original bhspec and kerrbb models. The 
lower abundance means that the absorption edge 
structure below 1 keV is less significant (~ 3% 
peak-to-peak), which makes the fit slightly better 
ixl = 1152/706 for M = 7 Mq and i = 67°), 
though the residuals are still clearly dominated by 
the edge structure below 1 keV. We also show the 
fit results based on the bhspec with 1/10-solar 
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abundance in table [Sj Fig. [3 and Fig. [5l Even 
though this is an underestimate of the metalhcity, 
the predicted edge features are still larger than 
those seen in the data (Fig. [7]), and thus the fits 
are worse than the KERRBB. 

This illustrates the key issue. The data in- 
dicate that the spectrum is broader than a sin- 
gle color temperature corrected, relativistically 
smeared model (kerrbb). Including ion opaci- 
ties as in BHSPEC gives broader spectra, but with 
the inevitable consequence edge features are pre- 
dicted. Since there is no evidence for these edges in 
the spectrum, incorporating sub-solar abundances 
alone in bhspec cannot explain why the model 
does not reproduce the data. 

Instead, this could point to the importance 
of including more atomic physics in bhspec. 
Currently it only uses bound-free opacity, but 
bound-bound (line) transitions could also be im- 
portant, especially if there is large scale turbu- 
lence/convective motion in the disk atmosphere. 
However, if these additional transitions primarily 
show up as absorption features, they will further 
enhance the depth of the atomic absorption, in- 
creasing the discrepancies. Instead, line and/or 
radiative recombination continuum emission could 
fill in some of the absorption features, while still 
retaining the broader continuum, though this 
seems somewhat fine-tuned. Self-irradiation of 
the disk may be a better mechanism, as this would 
drive the photosphere towards isothermality, re- 
moving the atomic features which do not appear 
to be present in the data but keeping the broader 
continuum. 

Whatever the origin of the fit residuals, their 
existence demonstrates the possibility of using the 
observational data to discriminate between diff'er- 
ent models of the vertical structure of the accre- 
tion disk, and motivates development of more so- 
phisticated disk atmosphere calculations. 

We are grateful to all the Suzaku team mem- 
bers. The present work is supported by Grant- 
in- Aid No. 19740113 from Ministry of Education, 
Culture, Sports, Science and Technology of Japan. 
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grants NSF AST-0807432, NASA NNX08AH24G, 
and NSF AST-0807444. 
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0.1 1 10 

Energy [keV] 

Fig. 5. — Comparison of the best fit bhspec 
and KERRBB (black dashed hne) with i — 67° 
and M — 7 Mq. The bhspec models are shown 
with 1-solar abundance (red solid line), 1/3-solar 
abundance (green dot-dashed line) and 1/10-solar 
abundance (blue dot line). Ratios of the best fit 
BHSPEC models to the KERRBB model are shown 
in the bottom panel. For the plotting purpose, the 
interstellar absorption and SIMPL are excluded in 
both models. 
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Energy (kcV) 

Fig. 6. — Same as Fig. 0] but for 1/3-solar abun- 
dance. 




Energy (kcV) 

Fig. 7. — Same as Fig.|4]but for 1/10-solar abun- 
dance. 
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